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ABSTRACT 

Emission from Active Galactic Nuclei is known to vary strongly over time over a 
wide energy band, but the origin of the variability and especially of the inter-band 
correlations is still not well established. Here we present the results of our X-ray and 
optical monitoring campaign of the quasar MR 2251-178, covering a period of 2.5 
years. The X-ray 2-10 keV flux is remarkably well correlated with the optical B, V 
and R bands, their fluctuations are almost simultaneous with a delay consistent with 
days and not larger than 4 days in either direction. The amplitude of variations 
shows an intriguing behaviour: rapid, large amplitude fluctuations over tens of days 
in the X-rays have only small counterparts in the optical bands, while the long-term 
trends over hundreds of days are stronger in the B band than in X-rays. We show that 
simple reprocessing models, where all the optical variability arises from the variable 
X-ray heating, cannot simultaneously explain the discrepant variability amplitudes 
on different time-scales and the short delays between X-ray and optical bands. We 
interpret the variability and correlations, in the optically-thick accretion disc plus 
corona scenario, as the result of intrinsic accretion rate variations modulating both 
X-ray and optical emission, together with reprocessing of X-rays by the accretion disc. 
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^ 1 INTRODUCTION 



Optical continuum variability is a universal property of 
radio-quiet AGN which has been studied for several decades, 
but its origin is still unclear. The optical continuum emis- 
sion almost certainly o riginates from the accretion disc 
l|Koratkar fc Blaed 1 19991 1 , so it is natural to assume that 
disc variability (e.g. through accretion rate fluctuations) 
drives the optical variability. Models where the variability 
is intrinsic to the disc, however, suffer from several prob- 
lems. Firstly, optical variability is seen on time-scales as 
short as a day, whereas accretion variability time-s cales in 
standard accretion discs (|Shakura fc Svunvaevlll973l ) should 
be long in the optically-emitting regions, comparable to 
the visco us time-scale which can be of the order of a year 
or more (jTreves et al.lll988r ). Secondly, variations are well- 
correlated in different optical bands, with only short (days) 
delays between bands, in the sense that longer-wavelength 
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varia t ions lag shorter-wavelength variations ([Wanders et al.l 
1 19971 : ICollier et aJ] I2001I : ICackett et al] |2007| ). Since local 
disc temperature T should decrease with radius R as T (x 
7?~^'*, shorter wavelengths originate from smaller radii. 
Therefore, in the case of inward- propagating fluctuations the 
long-wavelength variations lead, rather than lag, the short- 
wavelength variations. Also, the resulting delays would be 
large, comparable to the radial drift time-scale, contrary to 
the observed day-scale delays. 

To acc o unt for the observed optical properties, 
iKrolik et al.l l|l99ll ) suggested that the variability is driven 
by changes of the X-ray continuum which illuminates and 
thereby heats the disc, causing optical continuum variations. 
Since the X-ray emission is likely to be centrally concen- 
trated, blue bands should respond to the X-ray variations 
flrst, followed by red, with short delays corresponding to the 
light-travel time delay between the blue and red-emitting 
parts of the disc. The predicte d optical time-delay r scales 
with wavelength A as r oc A^''' (|Collier et al.lll999l ). and the 
observed dependencies of optical continuum lags on wave- 
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lengt h are consistent with this prediction IjCackett et alj 

A more stringent test of reprocessing models is to di- 
rectly compare X-ray and optical variations, using multi- 
wavelength monitoring campaigns. To date, a handful of 
these campaigns have been carried out, largely facilitated 
by the Rossi X-ray Timing Explorer (RXTE) in conjunc- 
tion with various ground based observatories. The re- 
sults have been mixed, with some AGN showing good 
evidence for correlated variability, as woul d be expected 
from reprocessing models (e.g. NGC 4 051, Peterson et al.l 
I2OOOI : IShemmeret"aIll2003l: N GC 5548, lUttlev et al.ll2003l : 
Mrk 509. lMarshall et al.ll200a) . while others show more com- 
plicated but p ossibly correlated behaviour (e.g. NGC 7469, 
iNandra et al.j 1998, 2000) and one AGN shows no apparent 
correlation at all (NGC 3516. lMaoz et al.ll2002l '). 

The best c ase for corre l ated v ariability to date is 
NGC 5548, but lUttlev et all (|2003l ) note that it is dif- 
ficult to reconcile the large fractional amplitude of opti- 
cal variability on time-scales of months-years with repro- 
cessing models, since the X-ray variability shows smaller- 
amplitude variations than the optical on these long time- 
scales. One would expect smaller relative variations in opti- 
cal since the intrinsic emission due to viscous disc heating 
would dilute any variable reprocessed component. Further- 
more, Gaskell (2007) notes that a simple energetics argu- 
ment rules out reprocessing as the dominant source of op- 
tical variability in many AGN (including NGC 5548) which 
have 'big blue bumps' dominating the total luminosity, sig- 
nificantly exceeding the X-ray luminosity available for re- 
processing. Detailed reproces sing calculations, produced by 
iKazanas fc NavakshinI l|200lf ) for the case of NGC 3516 show 
that simple reprocessing of X-rays on their own could not 
produce the observed optical variability. Also, in some AGN 
with simultaneous X-ray and optical monitoring, there is ev- 
idence that the optical may lead the X-ray variations, which 
might be expected if at least some of the optical variabil- 
ity is produced by int rinsic accretion fluctuations propagat- 
ing t hrough the disc l|Shemmer et al.l l2003l : [Marshall et al.l 
l2008l ). 

Clearly, no single model provides a satisfactory explana- 
tion for all the data. It is possible, however, that some com- 
bination of reprocessing and intrinsic accretion variations 
may explain the diverse range of optical/X-ray behaviour 
which is already observed in only a small sample of AGN. 
The location of the optical emitting region probably plays 
a key role in determining the balance of intrinsic versus re- 
processed variability. It is governed by the disc temperature, 
which scales inversely with radius R (in units of the Gravi- 
tational radius Rg = GM/(?) but also scales with black hole 
mass A/bh and accretion rate (as a fraction of the Edding- 
ton rate) rh as T oc [{m/MBii)R~^Y'^ JTreves et al.lll988l ). 
Thus over the 3 decade range in black hole mass expected for 
AGN, we expect the radius corresponding to a given tem- 
perature to change by a factor of 10 (at the same fractional 
accretion rate), perhaps even more for different accretion 
rates. Since disc temperature governs the radius w here peak 
optical emission is produced, lUttlev et al.l l|2003h suggested 
that a diverse range of optical/X-ray behaviour might result 
from the rang e of masses and accretion rates expected in 
AGN (see also lLi fc Caoll2008l '). For example, the AGN with 
the most massive black holes will have cool discs and very 



centrally concentrated optical-emitting regions, in terms of 
gravitational radii. Therefore, the disc variability time-scales 
will be short, comparable with the X-ray variability time- 
scales, so intrinsic disc variability as well as reprocessing may 
contribute significantly to the optical variations. In contrast, 
AGN with lower-mass black holes will have hotter discs so 
optical emission will originate from larger radii, where disc 
variability time-scales are very long compared to the corre- 
sponding time-scales in the innermost regions, so that only 
reprocessing may contribute significantly to the rapid vari- 
ability we observe. 

To test composite models of optical variability, it is 
necessary to carry out combined optical/X-ray monitoring 
of AGN with widely varying masses and accretion rates. 
To date, these campaigns have focused on Seyfert galaxies, 
which typically cover bolometric luminosities ranging from 
10'*'^ to lO**^ erg s~^ . Here we report the first optical/X-ray 
correlation for a more luminous radio-quiet AGN, the quasar 
MR 2251-178 which we have monitored simultaneously in 
X-ray and optical bands for the past 2.5 years. We describe 
the data in Section 2 and show the cross correlations in Sec- 
tion 3. In Section 4 we explore reprocessing and propagating 
fiuctuations scenarios to explain the X-ray/optical variabil- 
ity in this source and discuss the implications of our results 
in Section 5. 



2 DATA 

MR 2251-178 was monitored in the X-ray band using the 
Rossi X-ray Timing Explorer RXTE and in several opti- 
cal bands with the 1.3m SMARTS telescope in Chile and 
the Im telescope at the Wise Observatory. Below we give 
a brief description of the observational campaigns and the 
construction of the light curves. 



2.1 X-ray monitoring with RXTE 

We have monitored MR 2251-178 in the X-ray band with 
RXTE taking 1 ks exposure snapshot every 4.3 days. In the 
present paper, we include X-ray data from 2004 March 27 
to 2008 January 13. Within this period, an intensive moni- 
toring campaign was carried out with daily observations for 
three months, between 2007 August 2 and November 1. 

Data were obtained using the RXTE Proportional 
Counter Array (PC A), which is sensitive in the range ~ 
2.7 keV to 60 keV and consists of 5 Proportional Counter 
Units (PCUs). Since PCU has lost its xenon layer leading 
to a changed instrumental response and high background, 
and PCUs 1, 3 and 4 are regularly switched off for our obser- 
vations, we only extracted data from PCU 2. We use stan- 
dard good-time interval selection criteria (earth elevation, 
ELV> 10° and source pointing OFFSET< 0.02°). Back- 
ground data were created using the combined faint back- 
ground model and the SAA history file which are all current 
as of 2008 March. We extracted spectra in the 3-12 keV en- 
ergy range for each snapshot and used xspec to fit a simple 
absorbed power law model with the abs orption column fixed 
at th e Galactic value, A^'h = 2.8 x lO'^" (ILockman fc Savage! 
[1995I) to obtain an estimate of the 2-10 keV flux and its 
error. The top panel in Fig. [l] shows the resultant 2-10 keV 
flux light curve. 
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2.2 B and V band monitoring with SMARTS 

We used the ANDICAM instrument mounted on the 1.3m 
SMARTS telescope in Cerro Tololo, Chile, to take images ol 
MR 2251-178 in the B and V filters. The observational sam- 
pling was set up to match the RXTE monitoring, starting 
from August 2006. Observations were taken every four days 
between 2006 August 3 and 2007 December 14, with a daily 
sampling period in 2007 from August 2 to November 1. 

Two images were taken for each filter and observation 
date. We performed aperture photometry on MR 2251-178 
and two nearby, non-variable, reference stars in the field of 
view, of similar flux to the target, using an aperture radius 
of 10 pixels=3".7. The typical FWHM of the observations 
was l".5. We confirmed that the relative fiux of the refer- 
ence stars was constant throughout the campaign. Relative 
fiux light curves were constructed by dividing the fiux of the 
MR 2251-178 by the combined measurement of both stars. In 
all cases where both nightly observations were good, we com- 
bined the relative fluxes, to produce one light curve point 
per night of observation. A sample B band image is shown 
in Fig. [2] the target is marked by a circle in the Figure and 
both reference stars used are marked by squares and labelled 
SI and S2. The source appears point like in our images and 
the extraction radius is large enough that changes in see- 
ing produce a negligible effect in the measured flux. Using 
an aperture of 20 pixels produces an identical relative-flux 
light curve, with a flux offset of 3% compared to the 10 pixel 
aperture radius light curve. Visual inspection of optical spec- 
troscopic data obtained with the SMARTS 1.5m telescope 
(not shown here) does not reveal the presence of any absorp- 
tion features that could correspond to the underlying stellar 
population, so we conclude that the galaxy contibution to 
the AGN fiux is small. 



2.3 Calibration of the optical lightcurves 

The 1.3m SMARTS telescope observed Landolt standards 
fields every photometric night during the campaign to pro- 
duce zero-point magnitudes and extinction correction fac- 
tors for each filter. We used these zero-points and correction 
factors, together with the airmass of our observations, to 
calculate the magnitude of our reference stars. For the cali- 
bration, we performed aperture photometry on our reference 
stars with 20 pixel aperture radius, to match the aperture 
used on the Landolt standards, although the resulting mag- 
nitudes were almost identical to the values using 10 pixel 
apertures. The resulting average values, using 151 photo- 
metric measurements for B and 153 for V were B=14.00 
and V=13.71 for SI and B=14.90 and V=14.58 for S2, with 
errors (rms/number of points) of 0.003 mag. We converted 
these magnitudes into fluxes at 5500 A for the V band and 
4400A for B, assuming a flat spectrum (a = 0) within each 
band. We averaged the fluxes of both stars and used this 
value to cahbrate the flux hght curves of MR 2251-178. The 
error bars on all optical light curves correspond to the error 
on the relative flux and were calculated by propagating the 
error on the photometry of MR 2251-178 and the compar- 
ison stars. We do not include the error on the calibration 
zero-point of the reference stars because and error in this 
quantity produces a shift by the same factor in all the data 
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Figure 2. SMARTS field of view, the target MR 2251-178 and 
the two reference stars are labelled. The arrows pointing N and 
E are 1 arc minute long. 



points and does not contribute to the poin-to-point scatter 
in the light curves. 



2.4 B and R band monitoring with the Wise 
Observatory 

MR 2251-178 was monitored in the B and R optical bands 
between June 2005 and December 2007, with observations 
approximately every 2 weeks up to December 2006 and then 
weekly for the rest of the campaign. The data reduction and 
light curv e extraction were do ne following the procedure de- 
scribed in iKaspi et al.l ((20071) . As the fleld of view of these 
images is larger than in the SMARTS data, up to 11 non- 
variable reference stars could be used to create a relative flux 
light curve for MR 2251-178. We scaled the B band relative 
flux light curve to match the SMARTS B band data, over the 
periods of time where the light curves overlapped. To allow 
for differences in starlight contamination between the two 
light curves, possibly caused by the different apertures and 
average seeing, we fltted an additive offset as well as a rela- 
tive normalisation when matching the Wise Observatory to 
the SMARTS fluxes. The smaU best-fitting offset, 6 x 10"" 
erg s~^ cm~'^, was added to Wise Observatory light curve. 
For the R band calibration, we obtained the magnitudes of 
the reference stars SI and S2 published in the US Naval 
Observatory catalog (USNO VI) to correct the instrumen- 
tal magnitudes and to convert these into observed fluxes. 
The published R magnitudes are 12.82 for SI and 13.41 for 
S2. We used the calibrated flux of these reference stars to 
produce the light curve of MR 2251-178. 



3 X-RAY/OPTICAL CROSS CORRELATION 

As a first step to compare the variability in different energy 
bands, we cross correlated the light curves, using the discrete 
correlation function (DCF) method of lEdelson fc KroliM 
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Figure 1. MR 2251-178 light curves. From top to bottom: 2-10 keV X-rays; B band Wise Observatory data in blue open circles and 
B band SMARTS data in pink filled circles; V band SMARTS data, and R band Wise Observatory data. The X-ray light curve is only 
shown for dates with optical band coverage. Gaps in the light curves correspond to the epochs when the source is not visible from the 
ground. 



l| 19881 ) and the z-trans formed cross cor relation function 
(ZDCF) as described in lAlexande3 (|l997l '). The DCF mea- 
sures the degree of correlation between the two light curves 
as a function of time lag, i.e. displacement of one of the light 
curves on the time axis. A DCF value of 1 (-1) indicates com- 
pletely (anti)correlated data, a value of indicates the data 
are not correlated. Unless stated otherwise, a positive value 
of the time lag corresponds to X-rays leading. 

We first compute the ZDCF between the complete X- 
ray and B light curves as this method can deal with sparsely 
sampled data. The resulting correlation function has a cen- 
tral peak reaching a maximum correlation coefficient of 0.71 
at a lag of -1-9.6 days. This central peak is very broad, with 
correlation coefficient dropping to at lags of -142 and -|-212 
days. As is evident from the height of the ZDCF peak and 
from the light curves in Fig. [T] the long term behaviour of 
the X-ray and optical bands are well matched, but it is hard 
to obtain an accurate lag. To determine the correlation on 
shorter time-scales, we split the light curves into year-long 
segments and calculated the DCF independently for each 
segment, avoiding the yearly gaps. We later combined the 
resulting DCFs weighting each segment by the number of 
DCF points contained. The third segment of the light curves 
covers the intensive monitoring campaigns in B, V and X- 



ray bands. We re-sampled the intensive light curves taking 
one point every four days, to match the observational rate 
of the rest of the campaign. This was done to give similar 
weighting to each similar-length segment and avoid letting 
the intensive sampling period dominate the combined DCF. 

Figure |3] shows the DCF between the X-ray and B band 
light curves, calculated using the three segments shown in 
Fig- Hand including Wise Observatory and SMARTS B band 
data. The correlation has a broad central peak reaching a 
value of DCF=0.64 (DCF=0.72 for SMARTS data only) and 
a centroid lag of -f 3.6±9.3 days. The centroid was calculated 
as the weighted average of lags with corresponding DCF 
values greater than 80 percent of the maximum, within the 
central peak. We prefer the centroid to the DCF peak lag 
because it is less sensitive to the particular DCF binning 
used. The lags and errors were estimated using t he random 
sample selection method of lPeterson et al.l (1200J), selecting 
randomly 67 percent of the data points in the optical light 
curve together with their nearest X-ray data point and cal- 
culating the DCF centroid for 1000 such trials. The mean lag 
corresponds to the mean of the 1000 trial centroid-lags and 
the error corresponds to the rms spread of the trial centroid- 
lags. We repeated the calculation using B band SMARTS 
data only and obtained a virtually identical DCF. 
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Figure 3. The discrete correlation function (DCF) between the 
X-ray and B bands is plotted in black markers, positive lag val- 
ues correspond to X-rays leading. The dotted lines represent the 
95% extremes of the distribution of simulated X-ray light curves 
when correlated with the real B band light curve, the solid lines 
represent the 99% extremes of this distribution. The correlation 
peak around time lag=0 is significant above the 99% significance 
level. 
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Figure 4. DCF calculated using the daily sampled light curves 
only (MJD 54300 to 54400) between X-ray and B band in black 
dots and between B and V bands in red crosses. The X-ray vs 
V band DCF is very similar to the X-ray vs B band and is not 
shown for clarity. Positive lag values indicate higher-energy band 
leading. 



The X-ray light curve is equally well correlated with 
the V band. The DCF was calculated using the last two 
segments of the X-ray light curve and the corresponding 
segments of the SMARTS V band data, which was also re- 
sampled to match the sampling frequency of the rest of the 
light curve. The central peak reaches DCF=0.69 with a very 
similar shape and lag to the X-ray vs B band case. 

The DCF between X-ray and R band data was calcu- 
lated for the three light curve segments shown in Fig. [l] and 
later combined. The resulting X-ray vs R band DCF has a 
broad central peak reaching values of ~ 0.5 between lags 
of -24 and -1-80 days. The mean centroid lag is 29 days but 
the error on the lag is not well constrained. The significance 
of all the correlations was estimated using a Monte Carlo 
method described below in Sec. 13.11 



3.1 Testing the significance of the DCF 

We simulated X-ray light curves, uncorrelated with the ob- 
served data, and calculated the DCF between them and the 
observed optical light curve. By repeating this process many 
times, we estimate the probability of obtaining by chance 
equal or higher DCF p eaks than observed in th e real data. 
We used the method of lTimmer fc Konid ([19951) to produce 
red-noise light curves with a given underlying power spec- 
trum to generate the simulated X-rays. The model power 
spectrum was chosen to be a single-bend power law of slope 
equal to -1 at low frequencies, bending to a steeper slope 
at high frequencies, using parameters that match the shape 
and normalisation of MR 2251-178 X-ray power spectrum as 
measured by Summons et al. (in preparation). Light curve 
points were generated in bins of 0.1 days and for a length ten 
times longer than the observed X-ray light curve, in order 
to allow for the long-term trends which occur in red-noise 
data. The simulations were then sampled exactly as the real 
data and observational noise was added at the appropri- 
ate level. We calculated the DCF between these simulated 



light curves and the real optical data in the same way as 
was done with the real X-ray data. We repeated this pro- 
cess 2000 times, recording the DCF values. The mean, 95% 
and 99% extremes of the distribution of simulated B/X-ray 
band DCFs are plotted in dot-dashed, dotted and solid lines, 
respectively, in Fig. O As shown in the figure, the central 
peak in the real data DCF reaches higher than 99% of sim- 
ulated, uncorrelated light curves. The same is true for the 
X-ray vs V band DCF. In the case of X-ray vs R band, the 
DCF central peak reached the level of the top 95% of the 
simulations, but did not reach the 99% level. The R band 
light curve, however, is almost identical to the correspond- 
ing Wise Observatory B band light curve. Given that the 
total, Wise Observatory plus SMARTS, B band light curve 
is significantly correlated with the X-rays, it is likely that 
the drop in significance in the X-ray vs R correlation is only 
due to the smaller number of R data points available. 



3.2 Short time-scale variability 

We used the intensively sampled data taken between MJD 
54300 and 54400, to put tighter constraints on lags between 
the bands. The DCF between X-rays and B band and be- 
tween B and V bands, of the intensive sampling data only, 
are shown in Fig. |4l they both show a clear central peak 
around a lag of zero days. The error b ars in the Figure 
are th e standard DCF errors described in lEdelson fc KrolikI 
1 19881 ). The DCF centroids obtained for the main peaks are 
0.6 ± 3.1 days between X-ray and B, 1.3 ± 3.3 day for X- 
ray vs V and 0.5 ± 2.4 days between B and V. In all cases, 
positive lag values indicate higher-energy band leading. All 
the lags are consistent with days and are constrained to 
be less than ~ 4 days in either direction. The errors on the 
lags were again obtained using t he random sample selection 
method of lPeterson et al.l l|2004h . 
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Figure 5. X-ray, B and V light curves covering tiie intensive 
monitoring campaign. The B and V fluxes have been calculated 
assuming an effective bandwidth of 968 A for the B band and 
for the V band. 



3.3 Variability amplitudes 

Figure [S] shows the daily-sampled X-ray, B and V band light 
curves. The rapid X-ray fluctuations on time-scales of tens 
of days in the X-ray band have only small counterparts in 
the optical light curves, with a small or zero lag. As is clear 
from this figure and also Fig. [T] the long term trends over 
hundreds of days appear in all bands. The relative amplitude 
of these long-term fluctuations is larger in the B band than 
in the X-rays: the minimum to maximum flux ratio observed 
in the X-rays for the entire period shown in Fig.[l]is 2.8, be- 
tween the observations on MJD-50000=3915 to 4234, while 
the SMARTS B band data produce a ratio of 3.14 between 
days 3950 and 4240. An even larger flux ratio is shown on 
the period covered by the Wise Observatory data, when the 
B band flux drops by a factor of 4 between 3571 and 3950. 
Figure [6] shows the B band flux as a function of its near- 
est (less than 2 days apart) X-ray measurement, each nor- 
malised by the corresponding mean flux of the complete light 
curve. The dots correspond to the long term, '^ 2-week sam- 
pled campaign covering 900 days. Wise Observatory data; 
the squares represent a shorter, 4-day sampling light curve 
covering 500 days, and the crosses represent the intensive 
monitoring data, sampled daily for 90 days. The best fitting 
linear relation for each data set gets flatter for shorter time- 
scales probed: B/B = 1.07X/X-0.11 for the long term light 
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Figure 6. Flux-flux plot of the X-ray and B band data on differ- 
ent time-scales, for long (red dots), intermediate (blue squares) 
and short (black crosses) sampling patterns and light curve 
lengths. The best fitting linear relation for each data set gets 
flatter for shorter time-scales probed: B/B = 1.07X/X — 0.11 for 
the long term light curve (solid lines) B/B = 0.76X/X -f 0.24 for 
the intermediate (dashed line) and B/B = O.mX/X + 0.62 for 
the daily sampled, 90 day long light curve (dotted line). 



curve (solid hues) B/B = Q.7QX/X + 0.24 for the interme- 
diate (dashed line) and B/B = 0.39X/X -1-0.62 for the daily 
sampled, 90 day long light curve (dotted line). This change 
in slope shows how the B and X-ray light curves, although 
always correlated, have different time-scale behaviour, with 
rapid fluctuations being much stronger in the X-ray band 
and long term trends being stronger in the B band. 

The optical B and V band light curves from SMARTS 
are almost identical in shape and only differ in the amplitude 
of the fluctuations. Figure [7] shows the V flux as a function 
of B flux for the same epoch, for the evenly 4-day sampled 
data and for the intensive monitoring sample, separately. 
The flux-flux plots are almost linear, following a relation 

V = (0.398 ± 0.005)_B + (1.09 ± 0.03) x 10"^^ erg s"^ cm'^ 
(X^ = 137 for 79 dof) for the non-intensive segments and 

V = (0.375 ± 0.01)B + (1.31 ± 0.09) x lO'^^erg s'^ cm'^ 
{x^ = 93.4 for 73 dof) for the intensive sampling data. No- 
tice that the short term light curve has a slightly flatter 
flux-flux relation. Normalising each light curve to its mean 
before computing the flux-flux relation produces the relative 
changes of flux in each band: V/V = 0.75B/B + 0.027, i.e. 
the relative amplitude of V fluctuations is smaller than in B. 
This difference in relative amplitude can be intrinsic to the 
variability process but it can also be partly due to starlight 
contamination producing a stronger constant component in 
the V than B band light curves. 

The Wise Observatory B and R light curves are also well 
correlated following a linear relation oi R — 0.83B + 4.6 x 
10~^'^erg s~^ cm~^ shown in Fig. |S] Normalising each light 
curve by its mean produces a relative variability relation of 
R/R = 0.62B/B + 0.38. 
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Figure 7. Flux-flux plot of the SMARTS V and B band data, the 
crosses correspond to the 4-day sampling campaign and the open 
circles represent the intensive monitoring campaign. The flux in 
both bands varies almost proportionally, following a relation V = 
0.4S + 1.09 X 10~^^erg s"'^ cm~^, indicating that the B band 
varies proportionally more than the V band and that the latter 
contains a stronger constant offset. 
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Figure 8. Flux-flux plot of the Wise Observatory R and B band 
data. The flux in both bands varies almost proportionally, follow- 
ing a relation R = 0.83-B -I- 4.6 x 10~^^erg s~^ cm~^, indicating 
that the B band also varies proportionally more than the R band 
and that the latter contains a stronger constant offset. 



4 THE ORIGIN OF THE CORRELATION AND 
LAGS 

We have shown in the previous sections that the variability 
in the X-ray and optical B, V and R bands is well correlated 
and that any time lag between optical and X-rays is small 
(< 4 days). On time-scales of months, the fluctuations in the 
optical bands are similar to, and even larger than, the X-ray 
fluctuations while on shorter time-scales rapid, large ampli- 
tude X-ray flares have only small counterparts in the optical 
bands. In this section we will explore possible scenarios to 
explain the correlated variability between the bands. 



4.1 Reprocessing 

The small lags between optical and X-ray bands suggest 
their variations are connected through reprocessing of X- 
rays, possibly by an optically thick accretion disc, where the 
light travel time to the reprocessor is the main contributor to 
the time lags. This scenario could, in principle, also explain 
the lack of short time-scale fluctuations in the optical bands, 
as these can be smeared out by the finite light travel time 
to different parts of the reprocessor. 

To test the predictions of this scenario, we generated 
thermally-reprocessed light curves, using the observed X- 
ray light curves as inp u t. We used the method described in 
iKazanas fc NavakshinI (|2001f ) to construct the reprocessed 
light curves. This model assumes a source of X-rays on 
the axis of symmetry of the system, illuminating the ac- 
cretion disc from a height h in units of gravitational radii 
Rg = GM/(?. The optical light curve is constructed by 
adding the X-ray flux received at each location in the disc, 
which depends on the X-ray source height, to the locally dis- 
sipated gravitational energy, as a function of radius. Each 
annulus in the disc is then assumed to emit this total flux 
as a black body so the characteristic temperature and flux 
in a given optical band can be calculated. The light travel 
time to the different parts of the accretion disc and from 
there to the observer are taken into account to produce the 
expected delays between X-ray and optical bands. We left 
the disc accretion rate m as a free parameter, which governs 
the intrinsic, non-variable part of the optical flux. 

The observed X-ray flux was converted to total lumi- 
nosity assuming a distance of 271 Mpc to MR 2251-178, 
corresponding to its corrected redshift of z = 0.063, using 
Ho = 73 km/s/Mpc, Qm = 0.27 and nA=0.73. We followed 
the X-ray spectr al fitting of BeppoSAX data performed by 
lOrr et ai.l l|200ll ). to estimate the total (0.1-500 keV) X-ray 
flux incident on the disc and the reprocessed fraction. We 
considered a power law model with reflection off cold ma- 
terial with a small reflection fraction R. A fit with Xspec 
model pexriv, using i?=0.16 and ^ = 465.1 (|Orr et al.ll2"ooih 
produced a 0.01-500 keV fiux of 1.75 x 10"" erg s"^ cm'^, 
while the 2-10 keV flux was 4.64 x 10~" erg s"^ cm"^. 
Therefore, we multiplied the 2-10 keV light curve by a fac- 
tor of four to obtain the total incident X-rays. We also used 
the model fit to estimate the fraction of flux incident on a 
unit area of the thin disc which is absorbed, and therefore 
reprocessed into blackbody photons, obtaing a reprocessed 
fraction of 60 percent. 

The reprocessed optical light curves were scaled by the 
same distance of 271 Mpc to produce observed fluxes. We 
included an additional parameter A that multiplies the re- 
sulting optical flux to match the mean observed value. The 
function of this parameter is to absorb differences in pre- 
dicted and observed fluxes due to uncertainties in the dis- 
tance and extinction of the observed optical flux. 

We varied h, ifi, black hole mass, inner truncation ra- 
dius of the optically thick disc rin, inclination angle to the 
observer 9 and rescaling factor A to try to reproduce the 
observed B band light curve using the observed X-ray light 
curve as input. To first order, the effect of increasing h is 
to increase the lag between X-ray and reprocessed optical 
fluctuations and it also increases the disc covering fraction. 
The intrinsic disc accretion rate m governs the stable disc 
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flux, diluting the variability of the reprocessed component. 
Notice that the total optical luminosity depends on both m 
and the inner truncation radius of the disc. The inclination 
angle introduces smoothing of reprocessed fluctuations since 
the response observed on the farther side of the disc lags the 
response of the nearer side by a maximum time of twice the 
light crossing time of the truncation radius, the maximum 
being produced when the disc is viewed edge-on. The black 
hole mass deflnes the size of the gravitational radius, so all 
the distances and therefore time-scales scale linearly with 
this parameter. 

The model light curves corresponding to the best-fitting 
combinations of the parameters are shown in Fig.|9l together 
with the data for comparison. We first fitted each SMARTS 
B band light curve segment separately. The best-fitting pa- 
rameters for the first segment shown in Fig. [9] are h = 
1007^3, M = 10^MQ,m = 0.86, nn = 38i?5,A = 0.15 and 
for the second segment the parameters are h = 68-Rg , M = 
4.5 X 10*MQ,m = 0.38, nn = ?,ORg,A = 0.33, both with 
an edge-on view, 9 — 0.0. The high position of the X-ray 
source seems incompatible with an evaporating-disc corona 
but might correspond to the base of a weak jet. We ap- 
phed each set of parameters to the other segment to pro- 
duce the two pairs of reprocessed light curves, shown by 
the solid and dashed lines in the top panel of the figure. 
Clearly, the same set of model parameters cannot simulta- 
neously reproduce the shapes of the two light curves seg- 
ments and their relative offsets, which are governed by long- 
term trends. Within a single segment however, the model 
can reproduce the short term behaviour moderately well. 
This result immediately suggests that, while the short-term 
fluctuations may be caused by reprocessing, the long-term 
variability should have a separate origin, e.g. as intrinsic 
accretion disc fluctuations. 

To highlight the difficulty faced by simple reprocess- 
ing models in explaining the long-term optical variability, 
in the bottom panel of the figure we show a joint fit to 
both segments, with resulting parameters h — 100Rg,Ai[ = 
10^ M0,m = 0.58, nn = 83Rg,A = 0.46 and 9 = 0.0. 
The joint fit demonstrates that, in order to reproduce the 
long term fiuctuations, the model necessarily produces too 
much rapid optical variability, even at an edge-on angle 
of observation which maximizes the smoothing. Additional 
smoothing is produced by having a larger reprocessor i.e. 
pushing the inner truncation radius rin outwards, but even 
then the model cannot smooth out the rapid fiuctuations 
enough without introducing a noticeable optical to X-ray 
lag. Adding more intrinsic, constant, optical fiux to reduce 
the short-term fluctuations does not solve this issue as this 
would also reduce the long term trends, which are at least as 
strong in the B band as in X-rays. Therefore, X-ray repro- 
cessing cannot be the sole cause of all the observed optical 
variability. 

4.2 Propagating-fluctuations plus reprocessing 

A possible origin for the long-term optical variability is fluc- 
tuations travelling through the accretion disc. In this sce- 
nario the X-ray and optical long term trends can both be 
modulated by the accretion rate. A potential problem with 
this scenario is that if the accretion rate fluctuations prop- 
agate through a geometrically thin disc, their propagation 




B 


(N 


CJ 






o 


r/1 




£? 


X 

00 


OJ 








X 




P 




C 


^ 


pq 


6 




X 




^ 



4000 4200 

date [MJD-50000] 



4400 



Figure 9. Reprocessing: The solid and dashed lines represent the 
reprocessed B band light curves, calculated using the observed X- 
ray light curve as input. The dots represent the observed B-band 
light curve. In the top panel the reprocessing model has been 
fitted to only one segment at the time; the solid line corresponds 
to the fit to the first segment and the dashed line to the fit to the 
second. The solid line in the bottom panel shows the best fit to 
the two segments together. 



time can be long and we would expect to see time lags be- 
tween the bands, with lower energies leading. This problem 
can be circumvented if part of the X-rays are reprocessed 
by the disc. The large and rapid X-ray fluctuations can then 
produce some optical variability where the lower energy fluc- 
tuations will lag, therefore canceling part of the propagation 
effect. 

We constructed a propagating-fluctuation model that 
includes reprocessing of the X-rays to test whether it can 
reproduce the different amplitudes of variabilit y and cor- 
rect t ime lags. We followed the prescription of iLvubarskiil 
1 19971 ) ■ who proposes accretion rate fluctuations generated 
over a wide range in radii in the accretion disc, on the local 
viscous time-scale. The fluctuations propagate inward with 
viscous velocity through the disc and flnally modulate the X- 
ray emission at the centre. As the local time-scales decrease 
with radius, the more centrally concentrated emission con- 
tains the shortest time-scale fluctuations, seen especially in 
the X-rays. The resulting light curves for this model have a 
1// power spectrum, bending to steeper slopes above a char- 
acteristic frequency, which is directly related to the shortest 
variability time-scales included, i.e., the viscous time-scale 
at the inner edge of the ac cretion flow. In a standard disc 
(|Shakura &: SvunvaevligTa ) , this characteristic time-scale is 
proportional to the mass of the black hole and is a function 
of the disc thickness (H/R) and viscosity (a) parameters. 

We incorporated reprocessing into the implementation 
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of this model described in lArevalo fc UttlevI l|2006f ). by 
adding an X-ray source above the disc and on the axis of 
symmetry. We assumed that the optical flux is emitted ther- 
mally by the disc with local emissivity proportional to the 
radial gravitational energy release. This emissivity profile is 
first modulated by the propagating accretion rate fluctua- 
tions and then also by the variable flux received from the 
X-ray source. We incorporate the travel time of the accre- 
tion rate fluctuations to the central X-ray emitting region 
and the light travel time of the X-rays to the different annuli 
in the accretion disc. For simplicity, only a face-on viewing 
angle was considered. 

We used the observed X-ray flux and a distance of 27f 
Mpc as above, to obtain the X-ray luminosity. The black hole 
mass of MR 225f-f 78 has not been measured but, using the 
lower limit on the power spectrum break time-scale of Sum- 
mons et al. (in prep.) and the rel ation between this qua ntity 
and black hole mass derived in iMcHardv et al.l (|2006h . we 
derived a minimum mass of several times IQ^ Mq. We fixed 
the mass of the black hole to 8x fO*M0 and varied disc thick- 
ness and/or a parameter to obtain the correct X-ray power 
spectrum from the simulated X-ray light curve, obtaining 
a value of {H/R)^a = 0.08 at the innermost radius. As, in 
the model, the time-scales grow with black hole mass and 
decrease with a, these parameters are largely degenerate so 
we only picked one set that fits the data but note that other 
pairs would work equally well. We then fixed these parame- 
ters and varied the height of the X-ray source and intrinsic 
mean disc accretion rate, to produce different realisations of 
the optical light curves. 

One possible geometry for the accretion disc/corona 
system is a geometrically thin accretion disc that thick- 
ens towards the centre to produce a geometrically-thick, 
optically-thin corona. As the thickening of the flow produces 
shorter time-scale fluctuations at the same radius, the accre- 
tion rate that flnally modulates the X-rays can have fluctu- 
ations on time-scales much shorter than the optical emis- 
sion from the truncated thin disc. We kept the thickness 
to radius ratio H / R constant down to a truncation radius 
of Tt = fO-Rt,, with a value of {H/R)^a = 0.04, this pro- 
duces viscous fluctuations up to time-scales of 300 days at 
the truncation radius. Inside rt, we allowed H/R to grow 
with decreasing radius as {H/R)^a — 0.04(rt/r)^* to reach 
the desired value of {H/R)^a = 0.08 at r = 6Rg producing 
fluctuations on time-scales down to 5 days. We allowed the 
X-ray source to be higher than the top of the thick disc, 
but assumed that it was modulated by the thick flow fluc- 
tuations. The X-ray source height, however, had to be kept 
small at 4-6 Rg to produce sufficiently little reprocessing to 
reproduce the small relative size of the rapid optical fluctu- 
ations. 

Figure [TO] shows 500 day long simulated light curves. 
The top panel represents the fluctuations in the innermost 
region, corresponding to the X-ray light curve, the bottom 
panel shows a B band light curve modulated only by ac- 
cretion rate fluctuations and the middle panel shows this 
same B light curve when the effect of reprocessing of X-rays 
is incorporated. As expected, the resulting X-ray and opti- 
cal light curves are well correlated. Figure [TTI demonstrates 
the effects of reprocessing on the DCF, the dots show the 
DCF between the simulated X-ray and B band light curves 
if reprocessing is switched off (i.e. between the light curves 
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Figure 10. Simulated light curves of a propagating fluctuation 
model. The top panel represents the fluctuations in the innermost 
region, corresponding to the X-ray light curve, the bottom panel 
shows an accretion-disc emitted B-band light curve, modulated 
only by accretion rate fluctuations, the middle panel shows this 
same B band light curve when the effect of reprocessing of X-rays 
is incorporated. 



shown in the top and bottom panels in Fig.lfO|l. the B band 
leads the X-rays by approximately 50 days, which in this 
case is the travel time of the accretion fluctuations from 
the main B emitting region to the centre. If reprocessing is 
switched on, so that the rapid X-ray flares are imprinted 
as small optical fluctuations in the same light curve used 
above (shown in the middle panel of Fig.lfOp. the DCF peak 
shifts to a lag of zero days, represented by the crosses in 
Fig. If f I The exact shape of the DCF between B and X-ray 
light curves changes with different realisations of the sim- 
ulation even when exactly the same setup and parameters 
are used. In Fig. [T^we show a few examples of such DCFs, 
calculated for different light curve segments, each 500 days 
long. Given the changes in DCF peak height and direction 
of asymmetry, we cannot draw any conclusions from these 
structures in either the simulated or observed DCFs. 

The amplitude of short term fluctuations is strongly 
reduced in the optical compared to the X-rays and higher 
energy optical bands display larger variability than lower 
energy ones. The flux-flux relation of the simulated B and 
V light curves is similar to the real data, following a linear 
relation of the form V = 0.46B -I- f .9 x 10"^^ erg s"^ cm"^. 
In the model, the temperature radial profile of the accretion 
disc produces an emitting region that is more extended in 
the V than in the B band, therefore diluting more strongly 
in the V band variability. This causes the lower variabil- 
ity amplitude in V, especially on short time-scales, without 
including any starlight contribution that might dilute the 
variability further. 

The auto and cross correlation functions between op- 
tical light curves can also be reproduced. Figure [T3] shows 
the similarity of the DCFs between B and V bands of the 
real data and the simulated light curves shown in Fig. 1101 
These functions are almost symmetrical and very similar to 
the ACF, so we only display the DCFs for clarity. The sim- 
ulations for this particular realisation reproduce the right 
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Figure 11. DCF between the simulated X-ray and B band light 
curves shown in Fig. 1101 including reprocessing of X-rays (crosses) 
and not including reprocessing (dots). The additional B band 
fluctuations produced by reprocessing of X-rays move the DCF 
peak to lag. 
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Figure 12. Sample of DCFs between different realisations of sim- 
ulated X-ray and B band light curves, with the same input pa- 
rameters and length as those shown in Fig. 1101 Here we plot the 
DCFs as continuous lines for clarity. The main DCF peak ap- 
pears consistently around lag but the height of the peak and 
the direction of its asymmetry depend on the realisation. 



time-scale dependence of the variability and produce a lag 
between B and V light curves consistent with 0. 

A feature which this simple implementation fails to re- 
produce is the larger amplitude of long term optical fluctu- 
ations compared to the X-rays. This effect can appear for 
example if the fluctuations are damped as they propagate 
inward or as they transfer from the optically thick flow to the 
corona. As we only attempted a phenomenological approach 
to this fluctuating accretion model, we can easily introduce 
a damping parameter that reduces the X-ray long term fluc- 
tuations but this would not produce additional information 
on the physical processes involved. 



U 
Q 



' ' 


A' 


. 11 1 11 11 1 




&'+,» 


■ 




/■ .\ 


• B vs V real data 




: /: 


^ + B vs V simulated 


in 

d 

o 


/■ 
- /■ 


TL •• 




■ .+ 


+ - 


in 


- ++ 


++ - 


o 
1 


-/ 


¥ 




-*li 


+ 






+.< 




f 


i_ 


. — 1 


1 , , , , 1 , , , , 1 , 


... 1 .... 1 



-100 -50 

time lag [days] 



50 



100 



Figure 13. DCF between B and V band light curves for the real 
data (dots) and propagating-fluctuation plus reprocessing simu- 
lated data (crosses). 



5 CONCLUSIONS 

We have carried out a simultaneous monitoring campaign in 
X-ray and optical B, V and R bands on the quasar MR 2251- 
178 over two and a half years. All bands are significantly 
variable and their fluctuations are well correlated. The cross 
correlation functions show peaks around a lag of zero days, 
signiflcant over the 99% confidence level when compared 
to uncorrelated simulated data. All the delays between the 
bands are consistent with days. We used daily sampled 
light curves in X-ray, B and V to constrain the lags between 
these bands, obtaining values of 0.6 ± 3.1 days between X- 
ray and B, 1.3 ± 3.3 day for X-ray vs V and 0.5 ± 2.4 days 
between B and V lags, where positive values indicate higher 
energy band leading. The lag between X-ray and R band 
was determined using the long but more sparsely sampled 
light curve, obtaining a value of -4.5±16.8 days. 

The long term trends, over hundreds of days, are well 
matched by all the bands observed and are stronger in the 
B band than in X-rays, while the other optical bands show 
smaller-amplitude trends. The optical light curves are not 
corrected by host galaxy contribution, however, so the de- 
creased amplitude of variability in lower energy bands might 
be partly due to a constant galactic contribution. On short 
time-scales, of tens of days, there are large X-ray flares which 
appear strongly attenuated in the optical bands. The inten- 
sively sampled light curves show that the rapid X-ray flares 
do have optical counterparts, but that these are very weak. 

Pure reprocessing of X-rays cannot produce both the 
short and long time-scale optical variability: if the long 
term optical fluctuations were produced by reprocessing, we 
would expect similarly large rapid optical fluctuations, un- 
less these could be smoothed by light travel effects. Any 
geometry of the reprocessor, however, cannot smooth out 
the fluctuations on time-scales of 20-50 days without intro- 
ducing time-lags of a similar length, which are not observed. 
A sim ilar conclusion was reached bv lKazanas fc NavakshinI 
(|200ll ) for the Seyfert galaxy NGC 3516, where the X-ray- 
UV variability could not be explained by simple reprocess- 
ing. We show that reprocessing of the observed X-rays can- 
not reproduce the optical variability in MR 2251-178, either 
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producing rapid fluctuations that were too large or long term 
trends that were too small. The short term fluctuations how- 
ever could be reproduced reasonably well with the model if 
we limited the light curve lengths to ~ 100 days, which 
reduces the contribution of long term trends to each light 
curve segment. Therefore, a moderate amount of reprocess- 
ing acting on an already long-term variable optical emission 
can reproduce the observed behaviour. 

A simple way to decouple the amplitude of short and 
long term fluctuations is to have two distinct processes pro- 
ducing the optical variability, e.g. accretion rate fluctuations 
plus reprocessing. The covering fraction of the disc, together 
with the ratio of intrinsic disc emission to X-ray heating de- 
termine the fraction of optical flux arising from reprocessed 
X-rays. Therefore, the size of the rapid optical fluctuations 
depends on these parameters and can be largely indepen- 
dent of the amplitude of the long term trends. To test this 
possibility, we generated simulated light curves and com- 
pared their statistical properties with the observed data. In 
the model, fluctuations propagate inward through the ac- 
cretion flow, modulating the optical emission as they travel 
through the thermally emitting region and finally modulate 
the X-ray emission, assumed to be located in the centre. This 
process leads to the long term correlated variability. As the 
fluctuations are produced at every radius in the flow, on a 
radially dependent time-scale, the centrally emitted X-rays 
contain short time-scale fluctuations which are originally ab- 
sent in the optical light curves. The rapid fluctuations, only 
present in the X-ray emitting region are then introduced 
into the optical light curves through thermal reprocessing. 

Lastly, although propagation times across the accretion 
disc can be long and produce a large lag between optical 
and X-ray fluctuations, reprocessing reduces the lag and can 
even revert it to a value of the order of the light crossing 
time to the reprocessor. For a black hole mass of W^Mq 
the light crossing time of one gravitational radius is 5000 
s, so even for a large source height and inner truncation 
radius of 50 Rg the crossing time would be 4 days, less for a 
less massive black hole, which is within the observed limits 
on the lag. We show using the simulated light curves, that 
small-amplitude rapid fluctuations in the optical band can 
shift the CCF peak towards a lag of zero days, even when the 
long term trends of the optical lead the X-rays by ~50 days. 
As inward propagation of the fluctuations and reprocessing 
have opposite effects on the time lag, the relative amount of 
optical/UV variability produced by each process determines 
the sign of the lag. This might be the case in other monitored 
AGN where the optical/UV can either lead or lag the X-rays. 

An interesting feature of our fits to the optical vari- 
ability is that we expect the disc to subtend a small solid 
angle with respect to the X-ray source, in order to reduce 
the effects of X-ray reprocessing and produce the correct 
amplitude of optical rapid fluctuations. This small solid an- 
gle is consistent with the small amount of reflection inferred 
from the broadband X-ray spectrum of MR 2251-178 ob- 
served by BeppoSAX (|Orr et al.ll200in . which shows only a 
weak iron line with an equivalent width of EW ~ 70 eV, 
and modest reflection continuum (reflected fraction < 0.4). 
It is therefore possible that MR 2251-178 represents a source 
in a state where the thin disc is slightly truncated and/or 
the X-ray emitting region is at a small height over the disc. 



Optical/X-ray studies of other AGN should provide further 
useful constraints on their disc/corona geometries. 
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